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BURNING  OF  GRAPHITE  IN  AN  AIR  FLOW  AT  HIGH  TEMPERATURES 
N.  A.  Anfimov 


In  the  overwhelming  majority  of  theoretical  investigations  the 
examination  of  various  problems  of  the  boundary  layer  is  conducted 
with  the  use  of  a  presentation  of  component  mixtures  of  gases  as 
binary.  In  work  [1]  it  was  shown  that  this  presentation  gives  good 
results  and  is  fully  warranted  with  the  calculation  of  heat  exchange 
and  friction  on  an  impermeable  surface  in  dissociated  air. 

However,  the  presentation  of  a  multi-component  mixture  of  gases 
as  a  binary  proves  to  be  unsuitable  in  the  presence  in  the  mixture 
of  three  or  more  components  with  different  properties,  and  also  for 
a  pseudo-binary  mixture  of  gases  in  the  presence  of  interdiffusion 
within  each  group  of  components. 

The  importance  of  the  condition  of  constancy  of  direction  of 
diffusion  for  the  latter  case  can  be  demonstrated  in  the  example  of 
a  boundary  layer  on  the  surface  of  graphite  streamlined  by  high-tem- 
perature  dissociated  air. 

The  problem  on  the  breakdown  of  graphite  in  a  flow  of  highs=  tem¬ 
perature  air  is  also  of  independent  interest  since  graohite  is  a  oro 


Rising  material  for  use  at  high  temperatures  [2], 


1.  Scheme  or  bt'eakdown  of  graphite  in  a  flow  of  dissociated  air. 
.With  the  streamlining  of  a  graphite  body  by  dissociated  air  two  pro¬ 
cesses  occur  simultaneously  which  are  connected  with  carry-off  of  the 
mass  of  material  -  sublimation  and  chemical  interaction  between  carbon 
and  and  drag.  Let  us  examine  some  features  of  these  processes. 

Sublimation .  With  normal  pressure,  graphite  changes  to  a  gas¬ 
like  state,  by-passing  the  liquid  phase.  Graphite  vapors  consist  of 
various  gas-like  carbon  compounds:  C,  C2,  C^,  Cii,  C^,  Cg  and  [3]. 

In  the  case  of  a  substance  whose  vapors  consist  of  a  single  product. 


the  rate  of  evaporation  is  determined  by  the  Knudsen-Langmuir  ratio 

r  -  B <*»-*>  n  -1  v 

Gt  VS^RjWt  (1,1) 

Here  pi*  -  pressure  of  saturated  vapors  (equilibrium  pressure). 


f(.  -  pressure  of  the  vapors  of  the  material  at  the  surface,  a  -  ev¬ 
aporation  coefficient,  P.  -  gas  constant,  -  vapors*  molecular  weight, 
Tw  -  surface  temperature. 

Here,  it  turns  out  that  in  the  absence  of  chemical  reactions  in 
the  process  of  carry-off  of  mass,  the  pressure  of  vapor  at  the  surface 
is  close  to  equilibrium  [*!]. 

In  the  case  of  a  substance  whose  vapors  consist  of  several  gas¬ 
like  products,  .lust  as  in  the  case  of  graphite,  the  full  evaporation 
rate  is  determined  as  the  sum  of  evaporation  rates  in  the  form  of 
individual  products.  For  graphite,  according  to  the  data  in  work  [3], 
a0=O,3;  aCo=0. 5-1.0;  ac^O.l,  and  for  the  remaining  compounds  of  car¬ 
bon  the  evaporation  coefficient  is  close  to  null.  According  to  the 
data  of  other  authors,  we  obtain  other  values  of  evaporation  coeffi¬ 
cient  [5]. 
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Chemical  interaction  of  graphite  with  dissociated  air.  At  high 
temperatures  there  can  occur  the  following  chemical  reactions  between 
carbon  and  the  components  of  air  (we  have  written  only  independent 

reactions):  _  _ _ _ 

C+  O*fcC0,  C+  N^fcCN 
C  +  20 »» CO*,  '2C  +  2N^CiNa 

There  are  data  in  literature  [6]  which  attest  to  the  fact  that 
reactions  between  graphite  and  oxygen  are  heterogeneous.  This  means 
that  carbon  is  combined  with  oxygen  of  air  directly  in  the  solid 
phase,  by-passing  the  sublimation  stage.  Chemical  reactions  between 
graphite  and  oxygen  occure  in  several  stages  [7].  First,  atoms  and 
molecules  of  oxygen  diffuse  through  the  boundary  layer  to  the  graphite 
surface.  After  this,  they  are  adsorbed  by  the  surface.  Then  there 
occurs  a  strictly  chemical  reaction  between  carbon  and  oxygen.  The 
following  stage  -  desorption  of  products  of  the  reaction  from  the 
surface  and,  finally  -  diffusion  of  these  products  throughvthe  boun¬ 
dary  layer. 

The  speed  of  the  reaction  is  wholly  determined  by  the  speed  of 
the  slowest  of  the  elementary  processes.  At  low  temperatures  of  the 
surface,  this  determining  process  is  the  strictly  chemical  reaction. 
There  is,  presently,  a  sufficient  amount  of  data  on  the  kinetics  of 
burning  of  carbon  at  temperatures  up  to  1390-1*300°  K  [6],  These  data 


can  also  be  used  for  conducting  quantitative  calculations,  as  was 
dene  in  work  [3],  Here,  we  must  not  forget  that  the  constants  of  the 
rate  of  oxidation  of  graphite  strongly  depend  on  the  porosity  and 
on  technological  factors. 

At  temperatures  of  the  surface  on  the  order  of  3000°  K,  the  de¬ 
terminant  becomes  the  rate  of  diffusion  in  the  boundary  layer.  In  this 


case  the  chemical  composition  of  ras  in  the  surface  can  be  considered 
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equilibrium  since  the  gas-like  products  of  carry-off  continuously 
are  led  away  from  the  surface  by  means  of  convection  and  diffusion, 
through  the  boundary  layer.  Here,  the  rate  of  this  lead-off  depends 
on  the  interaction  of  various  mechanisms  of  the  carry-off  of  oass- 
ehenical  reactions  and  sublimation.  It  also  turned  cut  above  that 
reactions  between  carbon  and  oxygen  occur  on  the  surface  of  the  body; 
therefore,  it  is  incorrect  to  propose  that  graphite  first  sublimates 
then  enters  into  the  chemical  reaction  with  oxygen  as  was  done  in 
work  [9]»  In  this  case,  the  rate  of  the  process  is  limited  by  the 
rate  of  sublimation  of  carbon,  which  leads  to  a  decrease  in  the  rate 
of  carry-off  of  mass  and  an  Increase  in  temperature  of  the  decompo¬ 
sing  surface.  For  example,  according  to  the  data  in  work  [9],  for  the 
breakdown  of  spherical  body  made  of  graphite  with  a  diameter  of  2  m, 
with  streamlining  by  an  oxygen  flow  with  temperature  Ta=666 1°  K  and 
pressure  p=0.93  atm,  the  rate  of  carry-off  of  mass  is  acquired  as 
equal  to  3^0.013  kg/m"sec,  and  the  temperature  of  the  surface  ?v;% 

%3600°  H,  while  according  to  the  method  of  the  present  work  it  turns 
out  that  3=0.050  kg/m^sec  and  ?W=2515°K. 

hunever,  even  if  we  consider  the  reaction  of  carbon  with  oxygen 
tc  be  heterogeneous,  the  question  on  the  path  of  formation  of  gas-like 
carbon  remains  unclear  -  C,  C2,  0,,  Cj,,  C^,  Cg  and  Cy.  Carbon  vapors 
can  be  formed  as  a  result  of  sublimation  of  graphite  or  as  a  result 
of  dissociation  of  gas-like  products  of  the  chemical  interaction  of 
graphite  with  carbon  ana  nitrogen. 

We  propose  that  the  rate  of  dissociation  of  these  products  is  low 
in  comparison  with  the  rate  of  sublimation  and  that,  consequently,  the 
entire  gas-like  carbon  is  formed  as  a  result  of  sublimation  of  graphite. 


The  system  c  f  boundary  conditions  on  the  surface  of  graphite 


has  the  form: 

equation  of  continuity  for  carbon 


P  2V*(C>  +  Gv^^.Gy 

I  « 

equation  of  continuity  for  nitrogen 


p  2  CiVi  4"  Gy  2  Vi™  a  ==  0 

i  i 

equation  of  nonequilibrium  sublimation 

>  BCmIPc^—PcJ  '  » 

ecuation  of  equilibrium  for  indeoendent  chemical  reactions 


(1.2). 

(1.3) 

(l.H) 


;c?o. 


Po  = 

Ps  v 

Ps°  tr 

Vpo, 

Kv" 

Vpo.p».  * 

-r 

Pen 

=  Kr,.  PC*.  _  V 

P°a-K 

p?-Kr*+* 

Pc  Pa 

Fc'Prt* 

Pco 


=a:« 


(1.5) 


(1.6), 


u.y; 


curves  of  elasticity  of  saturated  carbon  vapors 

PC*H  =  h  (T »)  (*  =  1-4-7) 

balance  of  energy  on  the  surface  of  graphite 

—  ZpdVihi — Gy*.  +  Uvhlw  —  \t  ^ — mTJ  =  0 

In  the  equations  given  above,  p  -  density,  T  -  temperature,  h  - 
enthalpy,  X  -  coefficient  of  heat  conductivity,  -  partial  pressure 
of  the  i-components,  -  mass  concentration  of  the  i-component,  V*  - 
rate  of  diffusion,  -  mass  content  of  the  j-element  in  the  i-con- 

ponent,  G--  -  rate  of  carry-off  of  mass  of  the  material  in  the  gas¬ 
like  form,  ^  -  constants  of  reactions’  rates,  e  -  integral  radiating 
ability  of  graphite,  0  -  Stephan- Boltzmann  constant;  with  index  v:  we 
indicate  the  values  of  parameters  on  the  surface,  with  Index  1  -  the 
properties  of  solid  graphite. 

If  we  propose  that  graphite  breaks  down  exclusively  by  means  of 
sublimation  and  chemical  interaction,  then  from  the  equation  of  balance 
of  energy  we  can  determine  the  temperature  of  the  surface  of  graphite 


in  the  process  of  breakdown.  However,  witn  large  heat  flows  it  turns 
out  that  the  thusly-deternining  temperature  of  graphite  surface  can 
reach  too-high  values  (iJ000°K  and  more).  There  arise  doubts  relative 
to  the  ability  of  graphite  to  maintain  dynamic  forces  present  at  such 
a  temperature.  For  this,  in  order  to  acquire  the  higher  evaluation 
of  the  rate  of  breakdown  of  graphite  in  similar  conditions,  we  propose 
the  following  approach. 

We  propose  that  graphite  loses  mechanical  strength  with  some 
"threshold  temperature"  Ts.  Then,  if  from  equation  (1.7)  with  G*Gy 
it  follows  that  T...>TS,  then  we  consider  that  the  tenner  at.  ure  of  the 
surface  is  equal  to  its  maximum  value  T“,  and  the  excess  energy  occurs 
for  the  mechanical  breakdown  of  graphite,  in  which  regard  we  propose 
that  the  break-off  cf  solid  particles  of  graphite  occurs  without  ad¬ 
ditional  consumption  of  energy;  the  separated  particles  do  net  test 
the  physical-chemical  conversions  within  the  boundary  layer  and,  more- 


fusion  flows  which  approach  he  surface  of  the  body.  In  this  case, 
from  equation  of  balance  of  energy  (1.7)  we  can  determine  the  full 
speed  of  carry-off  of  mass  of  graphite  with  consideration  of  its  me- 
c hanical  breakdown . 

The  lead-off  of  heat  with  heat  conductivity  within  the  body, 
generally  speaking,  must  be  calculated  with  consideration  of  the  spe¬ 
cifics  of  the  oroblem  since  the  coefficient  of  heat  conductivity  of 


gr annate  is  extremeiy 


:.  We  will  orooose  that  the  heat  flow  withir 


graohite  corresnonds  to  the  stationary  linear  breakdown  of  the  semi¬ 


infinite  body.  In  this  case  (IhQ  -  enthalpy  of  unheated  graphite) 

TET  =  ®  (*»» — *w) 


and  equation  (1.7/  takes  the  form 


^  —  %f>CiV ift t — Gyhw — (G — Gy)  GAiji  —  tiTw*  =  0 


(1.8) 


After  the  system  of  boundary  conditions  on  the  surface  of  graph¬ 
ite  is  recorded,  for  a  full  solution  to  the  problem  we  must  have  a 
method  for  calculating  the  heat-  and  mass-exchange  on  the  surface  of 
the  breaking-down  body.  With  this  purpose,  we  can  use  either  appro¬ 
ximation  engineering  methods  or  methods  of  numerical  integration  of 
equations  of  boundary  layer. 

Whe  approximation  method  of  calculating  heat-  and  mass-exchange 
on  the  breaking-down  surface  leads  to  the  following:  a)  we  calculate 
the  coefficients  of  heat  exchange  on  an  equivalent  impermeable  sur¬ 
face;  b)  we  conduct  a  correction  for  the  injection  of  products  of 
breakdown  of  graphite  into  the  boundary  layer;  c)  we  calculate  the 
coefficients  of  mass  exchange  according  to  the  coefficients  of  heat 
exchange  with  the  use  of  an  analogy  between  the  processes  of  heat- 
and  mass  exchange. 

For  conducting  similar  calculations,  there  are  different  appro¬ 
ximation  formulas.  Some  of  them  are  done  in  works  [1G,11]  and  many 
others. 

A  numerical  solution  to  the  equations  of  multi- component  boundary 
layer  is  very  difficult  and  laborious,  although  it  permits  considering 
all  the  features  connected  with  the  complex  character  of  the  processes 
of  transfer  in  the  multi— component  boundary  layer  which  cannot  be 
fully  clarified  with  the  aid  of  approximation  solutions. 

In  the  present  work,  a  theoretical  investigation  of  the  carry-off 
of  mass  of  graphite  is  done  on  the  basis  of  a  numerical  solution  to  fch 
equations  of  a  multi-component  boundary  layer. 


2.  Numerical  solution  to  the  equations  of  a  mul ti -component 
boundary  layer  on  the  surface  of  graphite.  We  will  confine  ourselves 
to  an  examination  of  the  laminar  boundary  layer  in  the  area  of  the 
critical  point  of  an  ax i symmetrical  blunt  body.  In  this  case  the 
introduction  of  variables 


permits  switching  tc  normal  differential  equations.  In  these  variables 
the  system  of  equations  of  a  multi-component  boundary  layer  has  form 


vry +>f + 1(7- — /■*)  - '  (?,*) 

‘“If.  (2*23 

Here  x  and  y  -  axes  of  coordinates  directed  along  the  surface 
and  along  the  normals  to  it,  u0  -  speed  at  the  outer  boundary  of  the 
boundary  layer,  3  -  dimensionless  temperature,  vj.  -  dimensionless 
diffusion  flow  of  the  i-conponent,  c„  -  specific  heat  capacity  with 
constant  pressure,  u  -  viscosity,  f  -  dimensionless  function  of  fiov?, 

?  -  Prandtl  number,  w.:  -  speed  of  formation  of  the  i-component  as  a 
result  of  chemical  reactions.  With  index  e  v;e  indicate  the  values  of 
parameters  s,t  the  external  border  of  the  boundary  layer,  the  dash 
indicates  the  differentiating  with  respect  to  n. 

In  the  previous  work  of  the  author  [12],  devoted  to  the  carry-off 
of  mass  of  textolite  and  graphite,  the  diffusion  flows  were  calculated 
according  to  the  Wilke  approximation  formula  [13,  which  proposes  that 
the  diffusion  flow  of  each  component  is  determined  only  by  the  derived 
strict  concentration 


(2.U-) 


Here  D4  *  -  coefficients  of  binary  diffusion,  and  x*  -  mole  con- 
centration.  For  a  strict  formulation,  we  must  consider  the  dependence 


of  diffusion  flow  of  each  component  on  the  derived  concentrations  of  ! 
all  components  which  are  present  in  the  mixture.  The  relationships 
-between  derived  concentrations  and  diffusion  flows  are  presented  in 
work  [1*1];  after  transition  to  mass  concentrations  and  the  dimension¬ 
less  variables  accepted  in  the  present  work,  these  relationships  are 
written  in  the  form 

'  »  .  -  *  n  ~  —  -» 

•f  2  Ws  u=ii  +  7T  2  sr  ( 2  #  5 ) 

(J1J=|i/pDlj) 

rp 

Here,  -  Schmidt  number  for  binary  diffusion,  and  Dj_  -  di¬ 
mensionless  coefficients  of  thermodiffusion.  The  system  of  equations 

(2.5)  is  closed  by  relationship 

*!  +  •••  +  «»*  -»o.  (2,6) 

v/hich  shows  that  the  sum  of  diffusion  flows  is  identically  equal  to 
null  within  the  boundary  •  layer. 

A  direct  substitution  into  the  equations  of  continuity  of  the 
individual  components  (2.3)  of  diffusion  flows,  acquired  from  the  so¬ 
lution  of  the  linear  system  of  algebraic  equations  (2,5)  and  C2.6) 
leads  to  a  system  of  differentj  equations  of  the  second  order,  im- 
perrnissable  relative  to  the  higher  derivatives.  A  direct  solution  to 
this  system  with  iteration  methods  causes  objections  from  the  point 
of  tfiew  of  stability.  Therefore,  it  is  desirable  to  preliminarily 
solve  this  solution  of  equations  (examined  as  algebraic)  relative  to 
the  second  derived  concentrations.  On  the  whole,  this  method  proves 
to  be  ineffective. 

'lore  effective  is  another  method  which  is  used  below.  In  order 
to  obtain  a  system  of  differential  equations  for  concentrations,  we 
must  differentiate  the  system  for  diffusion  flows  [15,16],  Here,  the 
derived  diffusion  flows  are  expressed  with  the  aid  of  equations  of 
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continuity  of  individual  components,  and  the  second  derivative  of  ’ 
temperature  -  with  the  aid  of  energy  equation.  Here,  we  immediately 
acquire  a  system  of  differential  equations,  permissable  relative  to 
the  higher  derivatives  (with  an  accuracy  to  the  small  member  of  the 
set  (Mtci' i Mi)']  .  In  the  final  form,  the  equations  for  concentrations 
in  a  multi-component  boundary  layer  looks  like  this: 


c"  +  «<<?/  +  i><  «=  0,  Oi  =  — M  2'y- — fhi 

J-i  i 

,  N  c.V  N_1  w-i  .  i 

bi—c^MZjL)  Vi-/2  V/  +  2  +  (-2.7) 

i-i  ■  } 

+ ■?■{$[(?-)  +  V ~ S  **•<] 6 S  §7} + ( V) 


and  coefficients  K.y  and  are  determined  from  system  of  equations; 


w  ,, 
ei 


Ci'-cM S  g-  -  S  Vi- -Mine)' 


N-l 

2 

i-1  j-t 

With  a  practical  solution  to  the  system  of  equations  of  the  boun¬ 
dary  layer  in  the  present  work,  we  used  the  following  additional  pro¬ 
positions:  a)  the  rates  of  chemical  reactions  within  the  boundary 
layer  are  equal  to  null,  b)  the  contribution  of  thermodiffusion . to 
the  processes  of  transfer  of  mass  and  energy  can  be  disregarded. 

With  r  numerical  integration  of  the  full  system  of  equations  of 
a  multi-component  boundary  layer,  the  time  of  calculation  quickly 
grows  with  an  increase  in  the  number  of  components  (approximately 
proportionally  to  the  square  number  of  components).  Therefore,  we 
decided  to  stay  with  a  relatively  small  number  of  components,  with 
which  it  would  be  possible  to  conduct  a  rather  wide  series  of  calcu¬ 
lations  and  clarify  the  basic  features  of  carry-off  of  mass  of  graph¬ 
ite.  Thus,  the  number  of  components  in  the  boundary  layer  was  shortened 
to  eight:  0,  O2,  N,  1*2,  HO,  C,  CO  and  CM. 
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The  sequence  of  solution  to  the  system  of  equations  (2.1),  (2.2), 
(2,5)  -  (2.7)  with  consideration  of  the  boundary  conditions  (1.2)- 
(1,7)  was  as  follows, 

1.  As  the  null  approximation,  we  assign  random  values  of  the 
flow  function  on  the  wall,  temperatures  of  the  surface  and  concentra¬ 
tions  of  all  components  at  the  wall}  in  the  present  work,  the  null 
approximation  for  the  conditions  on  the  surface  was  selected  in  the 


following  manner: 


/»  -  o,  ‘  TV  -1000  +  0.25r, 
cco„  =  cn*,  =  0;5,  cqw  =  c0t¥l  =  eNo„  =  =  ccw  —  «cn*  —  0 


(2.8) 


2.  After  this,  we  assign  initial  profiles  of  speed,  temperature 
and  concentrations  within  the  boundary  layer.  For  this  purpose,  we 
used  profiles  similar  to  the  profile  of  speed  in  uncompressed  liquid 
(Blazius  profile). 

3.  Profiles  of  concentrations  were  differentiated,  and,  according 
to  values  c± ’  from  system  (2.5),  (2.6)  diffusion  flows  were  found. 

4 .  Then  system  (2.1),  (2.2)  and  (2.7)  was  solved  with  normal 
boundary  conditions  (2.8)  on  the  wall.  For  solving  this  system,  we 
used  a  method  of  trials,  the  use  of  which  in  equations  of  boundary 
layer  is  described  in  work  [17]. 

5.  According  to  the  acquired  solution  to  equations  of  boundary 
layer,  we  found  the  coefficients  of  mass  exchange  f?j_  and  coefficients 
of  heat  exchange  due  to  the  molecular  heat  conductivity 


v=V 


6.  On  the  basis  of  the  thusly-determined  coefficients  3^  and  a\ 
with  the  aid  of  boundary  conditions  on  the  breaking-down  surface  (1.2) 
-(1.7)  we  made  more  accurate  values  fw,  Tw  and  cj.  .  Here,  the  tie 
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between  value  of  flow  function  on  the  wall  and  the  speed  of  carry-off 


of  mass  of  the  material  in  gas  phase  is  given  by  expression 


If  it  turns  out  that  TW>TS,  then  thirs  is  an  indication  of  the 
presence  of  mechanical  carry-off.  In  this  case,  as  was  shown  above, 
we  propose  Tv;=T‘::‘ ,  and  from  relationship  (1.7)  we  find  the  full  speed 
of  carry-off  of  graphite  mass  with  consideration  of  its  mechanical 


breakdown. 


Table.  Parametex’s  of  potential  function  of  the  interaction  between 
particles  (asterisks  *  indicate  interpolation  according  to  molecular 


weights ) . 

KOMBCk 

—a 

M 

0 

2.96 

to 

44 

# 

0, 

3.433 

14 

113.2 

i14] 

N 

2.88 

*0 

38 

♦ 

N, 

3.681 

14 

91.46 

i**i  - 

NO 

3.47 

14 

119 

w  ■ 

C 

2.80 

33 

* 

CO 

3.59 

m 

110.3 

lM! 

CN 

3.50 

• 

100 

* 

Key:  1  -  component;  2  -  source. 

The  sequence  given  above  of  actions  was  repeated,  beginning  from 
point  3  up  to  the  time  when  the  relative  divergence  between  the  two 
sequential  approximations  did  not  become  less  than  some  value  5.  The 
approximations  were  compared  according  to  values  of  speed  of  carry- 
off  of  mass,  temperature  and  concentrations,  also  their  derivatives 
on  the  surface  of  the  body.  For  acquiring  convergence  with  an  accur¬ 
acy  of  S**10“i!  about  l|Q  approximations  were  required. 

The  method  of  calculating  the  properties  of  individual  components 
was  the  same  as  in  works  [1,  12,  15].  Here,  for  calculating  enthalpies 
of  individual  components  and  constants  of  equilibrium,  we  used  she 
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tables  from  work  [18].  It  was  proposed  that  the  interaction  of  all 
particles  present  in  the  mixture  is  described  by  the  Lennard-Johns 
potential  (6-12).  The  values  of  force  constants  used  are  given  in 
the  table. 

3.  Simplified  methods  of  numerical  solution  to  the  system  of 
equations  of  a  snul  ti -component  boundary  layer.  For  diffusion  flows., 


as  usual,  we  accepted  the  following  presentation 

etVl  =  —  D{^L 


(3.1) 


In  this  case,  the  equations  of  continuity  of  individual  compo¬ 


nents  are  simplified  and  have  the  form 


(v‘')+/c*,=o  (^pSy 


(3.2) 


where  S.;  -  Schmidt  number. 


The  diffusion  coefficients  in  (3.1)  are  determined  by  two  methods: 
by  the  Wilke  formula  (2.4)  and  by  the  formula  for  binary  mixture  of 
gases  where  the  coefficient  D_.  is  one  and  the  same  for  all  components. 

The  presentation  of  diffusion  coefficient  of  each  component  ac¬ 
cording  to  the  Wilke  formula  was  used  in  the  author’s  preceding  work 
[12],  In  this  case,  the  sum  of  diffusion  flows  dees  not  change  iden¬ 
tically  to  zero;  therefore,  with  the  solution  to  equations,  we  intro¬ 
duced  into  the  diffusion  flows  a  proportional  correction. 

In  the  case  of  the  presentation  of  a  multi- component  mixture  o,f 
gases  as  the  binary,  the  coefficient  of  diffusion  in  expression  (3.1) 
is  identical  for  all  components.  In  the  given  case,  the  mixture  of 
gases  can  be  divided  into  two  groups  with  approximately  identical 
properties:  1st  group  -  0,  13,  C;  2nd  group  -  Og,  CO,  CH. 

We  normally  recommend  as  the  single  coefficient  of  diffusion  taking 
the  coefficient  of  binary  diffusion  one  of  of  the  comoonents  of  the 


13 


iiwiiiiiiiiliiiiBiiniiii 


1st  group  with  one  of  the  components  of  the  2nd  group.  However,  we 
can  only  do  this  in  the  case  when  all  the  components  of  each  group 
are  differentiated  in  one  direction.  In  the  problem  of  carry-off  of 
graphite  mass,  the  latter  condition  is  not  accomplished.  So,  in  the 
1st  group  the  atoms  of  oxygen  and  nitrogen  diffuse  to  the  surface, 
and  the  atoms  of  carbon  -  from  itj  in  the  2nd  group  the  molecular 
oxygen  and  some  other  components  diffuse  to  the  surface,  and  the 
carbon  oxide  and  cyane  -  to  the  opposite  direction.  Therefore,  we  de¬ 
cided  to  conduct  calculations  for  two  coefficients  of  binary  diffu¬ 
sion:  <*)  O|  =  DO-C0.  <W  Di=°(VCO 

The  system  of  boundary  conditions  on  the  surface  of  graphite  and 
the  sequence  of  computations  with  the  presentation  of  diffusion  flows 
in  the  form  (3.1)  remain  unchanged. 

Figure  1. 


rigure  2. 


4.  Results  of  numerical  solution.  Calculations  were  conducted  .’o 
drag  with  parameters  Te=3030-S000oK,  p=5  atm. 

It  was  accepted  that  graphite  has  the  form  of  a  cylinder  with 
a  flat  end  with  a  diameter  of  Ik  mm.  The  integral  radiating  ability 
of  graonite  was  considered  eoual  to  0.3  [2],  evaooration  factor  -  a=l 


Numeric £.1  integration  of  equations  of  the  boundary  layer  was 
done  in  the  region  0<n<10  with  a  step  of  n  equau.  to  h=0.1.  The  control 
checks  ’with  step  h=0.05  showed  that  the  divergence  of  results  has  the 
same  order  as  the  assigned  degree  of  convergence  6. 

The  results  of  calculation  of  speed  of  carry-off  of  mass  G(kg/ 
/m^'sec)  and  temperature  of  the  surface  of  graphite  Tw  °K  depending 
on  temperature  of  braking  Te  °K  are  presented  in  Fig.  1  and  2,  where 
curve  1  -  accurate  solution;  curve  2  -  Wilke  formula  with  a  modifica¬ 
tion;  curves  3  -  binary  solution  for  D^=Dq  qq;  curves  4  -  binary  so¬ 
lution  for  cq.  bet  us  note  that  none  of  the  limitations  on 

value  of  temperature  of  the  surface  with  those  computed  apply. 

We  can  conclude  that  the  rate  of  breakdown  of  graphite,  which  is 
directly  connected  with  the  rates  of  diffusion  of  various  components 
within  the  boundary  layer ,  is  more  sensitive  to  the  method  cf  calcula¬ 
tion  of  diffusion  in  a  multi- component  boundary  layer.  The  temperature 
of  the  surface  is  more  stable  a  factor.  As  concerns  a  quantitative 
comparison  of  the  approximation  methods  with  an  accurate  solution,’* 
the  best  results  are  acquired  for  the  modified  Wilke  method.  With 
the  aid  of  this  method,  a  satisfactory  comparison  with  the  accurate 
colution  is  obtained  for  the  entire  studied  range  of  change  in  para¬ 
meters,  in  which  regard  the  error  of  determining  the  speed  of  break¬ 
down  with  Te=30G0-40Q0°K  is  less  than  l£,  and  with  Te=500Q-6000°K  the 
error  increases  to  53. 

The  fact  that  the  Wilke  approximation  formula  (with  subsequent 
modification)  qualitatively  accurately  describes  the  process  of.  dif¬ 
fusion  in  a  multi-component  boundary  layer  on  breaking-down  graphite, 
we  can  see  from  a  comparison  of  values  of  the  coefficients  of  mass- 
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exchange  3^,  computed  for  various  components  with  p=5  atm  and  T^= 
=5000°K  conforms  to  the  accurate  solution  (3i=3i’)  and  according  to 
the  Wilke  formula  with  modification  (3i=3^" ) 

c  CO  CM:  O  0.  K  N,  HO  ' 

0,'=  1.812  1.208  1.253  1.508  1.M5  1.591  0.973  1.145 
0/  =  1.856  1.259  1.297  1.418  1.105  ^1.533  1.150  1.106 

For  all  components,  with  the  exception  of  molecular  nitrogen, 

the  divergence  does  not  exceed  G%.  In  molecular  nitrogen  the  diver¬ 
gence  reaches  133,  in  which  regard  this  circumstance  is  connected 
with  the  anomolous  behavior  of  the  coefficient  of  mass  exchange  for 
■I which  can  be  established  only  as  a  result  of  solving  the  full 
system  of  differential  equations  and  is  not  transferred  by  approxi¬ 
mation  methods  (binary  mixture,  Wilke  formula).  With  a  calculation 
of  the  boundary  layer  with  the  use  of  presentation  of  binary  diffu-^ 
sion,  the  coefficients  of  mass  exchange  were  identical  for  all  the 
components  and  withe  the  indicated  conditions  composed,  respectively, 
3~1500  and  S^l^TS  for  two  values  (which  we  used)  of  coefficient  of 
binary  diffusion  Dq  qq  and  D q2  Cq. 

The  use  of  the  presentation  of  a  multi- component  mixture  of  gases 
as  a  binary  does  not  permit  the  satisfactorily-accurate  consideration 
of  all  the  laws  of  change  in  speed  of  breakdown  of  graphite  with  a 
change  in  temperature  of  impact.  However,  ‘with  the  aid  of  a  properly 
selected  coefficient  of  binary  diffusion,  we  can  obtain  a  satisfactory 
coincidence  at  one  or  another  temperature.  For  example,  with  Te=3000°K, 
the  difference  from  the  exact  solution  is  1.3$  for  bj_=Do2  qq»  and  with 
Te=6000°K,  the  difference  from  the  exact  solution  is  2!?  for  Dj_=Dj  C0» 

If  an  identical  coefficient  of  binary  diffusion  is  used  in  the  entire 


range  of  change  in  impact  temperature,  the  error  in  calculating “the 
rate  of  breakdown  of  graphite  can  reach  253.  We., can  conceive  an  ap- 
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proximation  method  of  calculating  a  multi-component  boundary  layer 
based  on  the  introduction  of  a  "determining  coefficient  of  diffusion" 
similar  to  the  method  of  the  determining  temperature  or  determining 
enthalpy,  with  the  aid  of  which  the  results  for  heat  exchange  or  fric¬ 
tion  for  uncompressed  liquid  with  constant  properties  are  transferred 
to  the  case  of  compressed  gas. 

As  we  noted  above,  the  temperature  of  the  surface  of  graphite 
depends  on  the  method  of  calculation  of  the  rates  of  diffusion  less 
strongly.  So,  in  the  investigated  range  of  parameter  change,  the  max¬ 
imum  divergence  in  value  of  surface  temperature  does  not  exceed  100°K, 
i.e.  3%. 


An  increase  in  rate  of  carry-off  of  mass  of  graphite  with  an 
increase  in  impact  temperature  of  the  air  flow  is  explained  by  two 
factors.  First,  we  see  here  an  increase  in  the  degree  of  dissociation 
of  oxygen,  as  a  result  of  which  the  diffusion  flow  of  oxygen  to  the 
surface  grows ,  since  the  atoms  of  oxygen  possess  larger  diffusion 
coefficients  than  molecules.  Secondly,  there  occurs  here  a  change  in 
the  mechanism  of  carry-off  of  mass  of  graphite.  Since  tfith  an  increase 
in  impact  temperature  the  surface  temperature  grows  simultaneously, 
then  besides  the  formation  of  a  carbon  oxide  the  role  of  sublimation 
of  carbon  (at  low  pressures)  and  formation  of  cyanogen  (at  high  pres¬ 
sures)  begins  to  increase.  Figure  _  shows  the  relative  fraction  of 
sublimation  of  graphite  Vq=Oq/0  and  tne  formation  of  cyanogen  Vqjj3 
=:1CnC”/’-CNG  in  Process  °-  carry-off  of  mass  of  graphite.  VJe  can 
see  that  with  T  =6000°K  and  o=5  atm,  dur  to  the  formation  of  cyanogen 
more  than  30%  of  the  mass  of  graphite  is  carried  off. 

Figure  shows  the  influence  of  limitation  of  surface  temperature 
of  graphite  on  its  mechanical  breakdown.  A±ong  the  v-axis  lies  the 
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fraction  of  carry-off  of  graphite  mass  in  the  gas-like  form,  and 
along  the  x-axis  -  maximum  surface  temperature.  Calculations  were  Con¬ 
ducted  for  the  following  parameters  of  drag:  p=5  atm,  Te=60C)0°K. 

The  results  of  calculations  show  that  the  slight  decrease  in  surface 
temperature  of  graphite  as  a  result  of  loss  of  mechanical  strength 
leads  to  a  sharp  increase  in  mechanical  carrv-off  of  mass.  We  can  also 
make  an  inverse  conclusion:  if  the  fraction  of  mechanical  breakdown 
of  graphite  is  not  very  great,  its  surface  temperature  lier^  hardly 
decreases  at  all. 

The  degree  of  the  nonequilibrium  state  of  graphite  sublimation 
for  the  mechanism  of  mass  carry-off  used  proved  to  be  very  low.  The 
pressure  of  carbon  vapors  can.  vrith  an  accuracy  sufficient  for  prac¬ 
tice,  be  considered  equal  to  the  pressure  of  saturated  vapors  since 
it  differs  from  this  maximum  value  by  only  1-2?. 


Figure  3.  ft\ 
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